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We study anisotropic lattice strips in the presence of a magnetic field in the quantum Hall effect 
regime. At specific magnetic fields, causing resonant Umklapp scattering, the system is gapped in 
the bulk and supports chiral edge states in close analogy to topological insulators. In electron gases 
with stripes, these gaps result in plateaus for the Hall conductivity exactly at the known fillings n/m 
(both positive integers and m odd) for the integer and fractional quantum Hall effect. For double 
strips we find topological phase transitions with phases that support midgap edge states with flat 
dispersion. The topological effects predicted here could be tested directly in optical lattices. 
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Introduction. Condensed matter systems with topo- 
logical properties have attracted wide attention over 
the years.— ~— E.g., the integer and fractional quantum 
Hall effects (IQHE and FQHE^ find their origin in 
the topology of the system.— ~— Similarly, band insu- 
lators with topological properties have become of cen- 
tral interest recently^^^ as well as exotic topological 
states like fractionally charged fermionsiSr— or Majorana 
fermions <2&~— 

Here, we study two-dimensional (2D) strips in mag- 
netic fields, both analytically and numerically, modeled 
by an anisotropic tight-binding lattice. We identify a 
striking mechanism by which the magnetic field induces 
resonant Umklapp scattering (across Brillouin zones) 
that opens a gap in the bulk spectrum and results in chi- 
ral edge states in analogy to topological insulators. Quite 
remarkably, the resonant scattering occurs at well-known 
filling factors for the IQHE^ and FQHE& v = n/m, where 
n, m arc positive integers and to odd. We argue below 
that this mechanism could shed new light on the QHE for 
2D electron gases as well, where the formation e.g. of a 
periodic structure (energetically favored also by a Peierls 
transition) might support the periodic structure needed 
for the Umklapp scattering. 

Finally, we consider a double strip of spinless fermions, 
or, equivalcntly, a single strip with spinful fermions. 
Here, we find two topological phase transitions accompa- 
nied by a closing and reopening of the bulk gap, and, as a 
result, three distinct phases. The trivial phase is without 
edge states. The first topological phase is similar to the 
one discussed above and carries two propagating chiral 
modes at each edge for v = 1. The second topological 
phase has only one state at each edge. Quite remark- 
ably, its dispersion is flat throughout the Brillouin zone, 
making this phase an attractive playground for studying 
interaction effects. 

Anisotropic tight-binding model. Wc consider a 2D 
tight-binding model of a strip that is of width W in x- 
and extended in y-direction, see Fig. [TJt. The unit cell 
is composed of two lattice sites (a = ±1) along y that 
are distinguished by two hopping amplitudes, t y i and t y <z- 
Every site is labeled by three indices n,m, and a, where 
n (to) denotes the position of the unit cell along the x- 
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FIG. 1. (a) Strip: two-dimensional lattice of width in x- 
direction, W , with unit cell defined by the lattice constants 
a x and a y . The hopping amplitudes in y-direction, t v \ and 
t y 2, carry the phase <f>, arising from a perpendicular magnetic 
field B, and we assume t x 3> t y i,t y 2- (b) Doubly degenerate 
spectrum of H x [see Eqs. Q] for the rows with a = 1 (upper 
blue) and with a = 1 (lower green). The hoppings t y i (dashed 
line) and t y 2 (dotted line) induce resonant scattering between 
right (R a ) and left (L&) movers, which open gaps at the Fermi 
wavevectors ±/cf defined by the chemical potential fi. 



(y-) axis. The hopping along x is described by 

H x = —t x i C n+l,m.(7 C n,m,a + h.C.), (1) 

n,m,o 

where t x is the hopping amplitude in x-direction and 
c n .m,a the annihilation operator acting on a spinless 
fermion at site (n, to, a), and the sum runs over all sites. 
The hopping along y is described by 



+ ty^cl m+hl c n ^ ri + h.c.), (2) 

Without loss of generality, we consider t V 2 > t y \ > 
0. The phase <f> is generated by a uniform magnetic 
field B applied in perpendicular z-direction, see Fig. 
[T] We choose the corresponding vector potential A, to 
be along the y-axis, A = (Bx)e y , yielding the phase 
<f) = eBa x a y /2hc. Here, the corresponding lat- 

tice constants. 

Chiral edge states. Taking into account translational 
invariance of the system in y-direction, we introduce 
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the momentum k y via Fourier transformation, see Ap- 
pendix [A] The Hamiltonians become diagonal in k y , 

Le -> H * = -^T,n,k y A C l+hk y ,a C n,k y ,a + H.C.), and 

/i.e.]. Thus, the eigenfunctions of H = H x + H y fac- 
torize as e %hyV ^ v {x), where we focus now on ipk y (%) and 
treat k y as parameter. 

Assuming for the moment periodic boundary condi- 
tions also in x-direction, we introduce a momentum k x , 
see Appendix [AJ Immediately, the well-known spectrum 
of H x follows, 6 (j = — 2t x cos(k x a x ), which is twofold de- 
generate in a. The chemical potential /i is fixed such that 
the Fermi wavevector hp is connected to the phase by 
(j) = 2kpa x . Next, we allow for hopping along y as a small 
perturbation to the x-hopping, i.e., t x 3> t y \,t y 2, see Fig. 
[TJ3. To obtain analytical solutions, it is most convenient 
to go to the continuum descriptio n 26 ! 39 . The annihila- 
tion operator 'I'(x) close to the Fermi level can be repre- 
sented in terms of slowly varying right [R a {x)] and left 
[L c (x)\ movers, $(x) = ~£<r R<x{x)e ikpx + L a {x)e- lkFX . 
The corresponding Hamiltonian density H can be rewrit- 
ten in terms of the Pauli matrices Tj (<Xj), acting on 
the right-left mover (lattice) subspace (see Appendix iBj) 
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FIG. 2. Spectrum E(k y ) of left edge state (red line or dots) 
propagating along y for a strip (t y i/t x = 0.02, t y ijt x = 0.1) 
of width W/a x — 801 and with phase (f> = 7r/2, obtained (a) 
analytically [see Eq. Q] and (b) numerically [see Eqs. (fJJ) 
and p])]. For k~ < k y < k+, there exists one edge state 
at each edge. The left (red dots) and the right (blue dots) 
edge states are chiral and propagate in opposite ^-directions. 
For each k v [(c) k y a v = n, (d) k v a v = 137r/12] there is one 
left (red dots) and one right (blue dots) edge state if Eq. 
([5]) is satisfied. Here, e(N) corresponds to the Ath energy 
level. The probability density |i/v| 2 (e) [(f)] of the left [right] 
localized state decays exponentially in agreement with the 
analytical result, Eq. (|C1|) in Appendix [C] 



* = (Ri,Li,Ri,L{), as 

H = Hv F kT 3 + "Y^m + cr 2 T 2 ) + -y 



(crrn - cr 2 r 2 ) 



x cos(k y a y ) + (ct 2 ti + ctit 2 ) sm(k y a y ) 



(3) 



Here, fi-fc = —ihd x is the momentum operator with eigen- 
values fc taken from the corresponding Fermi points ±/cf, 
and vf = 2(t x /H)a x sm(kpa x ) is the Fermi velocity. The 
spectrum with periodic boundary conditions in x- and 

y-directions is given by ej,± = ±4/ (Hvpk) 2 + t 2 ^ where 

I = 1,2. This mechanism of opening a gap by oscilla- 
tory terms causing resonant scattering between the Fermi 
points is similar to a Peierls transition^ Next, we turn 
to a strip of finite width W, see Fig. [TJ We note that the 
bulk spectrum ej ± is fully gapped, so states localized at 
the edges can potentially exist. To explore this possibility 
we consider a semi-infinite nanowire (x > 0) and follow 
the method developed in Refs. [26| and Unassuming that 
the localization length of bound states £ is much smaller 
than W. This allows us to impose vanishing boundary 
conditions only at x = 0, ipk y {x)\x=o = (i>i,i>i)\x=o = 0. 
This boundary condition is fulfilled only at one energy 
inside the gap \E\ <t y \, 



ty\t y 2 sm(k y a y ) 



E(ky) = —r 

+ *»2 - 2t„i£„ 2 cos(k y a y ) 

if the following condition is satisfied, 

tyi > t y 2 cos(k y a y ). 



(4) 



(5) 



The edge states exist for momenta k y E (fc_, fc + ), where 
k±a y = 7r ± arcsin(i y i/i y2 ). An edge state touches 
a boundary of the gap at k± and afterwards disap- 
pears in the bulk spectrum of the dclocalizcd states, see 
Fig. [2] The only regime in which the edge state ex- 
hibits all momenta corresponds to the uniform strip with 
t y i = t y 2- The localization length £ is determined by 

£ = hvp / ^Jt 2 yl — E 2 , with wavefunction given in Ap- 
pendix [Cj The edge state gets delocalized if its energy 
is close to the boundary of the gap, so that £ becomes 
comparable to W . Similarly, we can search for the so- 
lution decaying to the right, x < 0, and obtain Eq. ^ 
with reversed sign, E(k y ) — > —E(k y ). 

We have confirmed above results by diagonalizing the 
tight-binding Hamiltonian H (in fc^-representation) nu- 
merically, see Fig. [2j The spectrum E(k y ) of the edge 
states localized along x and propagating along y shows 
that at any fixed energy inside the gap there can be only 
one edge state at a given edge, see Fig. [3J Moreover, 
the edge states are chiral, as can be seen from the ve- 
locity, v = dE/dk y , which is negative (positive) for the 
left (right) edge state. This means that transport along a 
given edge of the strip can occur only in one direction de- 
termined by the direction of the B-field, see Fig. [TJ Quite 
remarkably, the obtained spectrum of edge states is of the 
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same form as for topological insulators^i with a sin- 
gle Dirac cone consisting of two crossing non-degenerate 
subgap modes. Due to the macroscopic separation of op- 
posite edges, £ <C W, these modes are protected from 
getting scattered into each other by impurities, phonons 
or interaction effects, so that the Dirac cone cannot be 
eliminated by perturbations that are local and smaller 
than the gap. Thus, the edge states are topologically 
stable. 

Umklapp scattering. We note that the system consid- 
ered here is equivalent to a 2D system in the QHE regime. 
The above choice of magnetic field corresponds to the 
IQHE with filling factor v = 1, which is in agreement 
with one chiral mode at each edge. To explore the pos- 
sibility of inducing quantum Hall physics at other filling 
factors, we fix the chemical potential \x = —y2t x , so that 
the system has local particle-hole symmetry, and change 
the B-field. Above, the phase <j>, generated by the mag- 
netic field, was equal to 7r/2 for kp = tt /Aa x . However, 
this is not the only choice of phase leading to the opening 
of a gap A fl at the Fermi level. Due to the periodicity 
of the spectrum, resonant scattering between branches of 
e a occurs also via Umklapp scattering between different 
Brillouin zones, with a phase = + 2irq, where q 
is an integer, n a positive integer, and p a positive odd 
integer with p < 2n and coprime to n. As a result, the 
Fermi level lies in the bulk gap for the filling factors^ 
v = n/(4qn±p), which can be rewritten as v = n/m, 
where m > is an odd integer. The size of the gap can 
be estimated as A g oc t v i(t y i /t x )^ n ~ 1 ^ (assuming for sim- 
plicity t y x = t y 2). Finally, we remark that we checked 
numerically that the gap A g never closes for any finite 
ratio of t x and t v i larger or smaller than one. 

FQHE in 2D electron gas. We conjecture that the same 
mechanism of resonant Umklapp scattering can also lead 
to the integer or fractional QHE in 2D electron gases. 
At high magnetic fields interaction effects get strongly 
enhanced and electrons tend to order themselves into 
periodic structures ! 15 ! 43 " — In particular, we assume the 
formation of stripes that are aligned along x and peri- 
odically repeated in y. While particles can hop between 
stripes, they move now continuously inside them with 
quadratic dispersion. Thus, the perturbative solutions 
found above in terms of right- and left-movers still apply. 
In addition, we assume that the interaction generates a 
charge-density wave at wavevector K inside the stripe, 
providing an effective periodic potential in x, which will 
lead to gaps. Thus, K becomes the period of the Bril- 
louin zone, and, at 1/4-filling of the lowest subband, we 
have K = 8k F £L 

Again, the B- field leads to a gap at kp only if it re- 
sults in phases commensurable to kp, i.e. eBa y /2hc = 
±2kF^ + qK, which is equivalent to v = n/ (Aqn ± p) = 
n/m. In this regime, there is an additional energy gain 
due to a Peierls transition, favoring even more a forma- 
tion of periodic structures with gaps. Moreover, from 
this mean field scenario it follows that the IQHE is 
more stable against disorder than the FQHE since the 



latter requires Umklapp scattering through higher Bril- 
louin zones. The gap and the edge states can be tested 
in transport experiments. For example, the Hall con- 
ductance an exhibits plateaus on the classical depen- 
dence curve o~h oc l/B, if the Fermi level lies in the 
gap. This can be shown by using the Streda formula^ 
o~h = ec j where n is the bulk particle density 

which is uniquely determined by the magnetic field via 
the relation veBa y /2hc — 2kp (for this it is crucial that 
K depends on kp). If fj. lies in the gap opened by the 
Umklapp scattering, the change in the density for fixed 
(i, dn, due to a change in the magnetic field, dB, is 
given by dn — ^ F j 2 = v ^dB. Hence, the conduc- 
tance assumes the FQHE plateaus, oh = ve 2 /he, with 
v = n/m and independent of any lattice parameters. 
The width of the plateaus is determined by the gap size 
A g oc tyi{t y il p) i - n - 1 \ We note that the FQHE can be 
mapped back to the IQHE by redefining the charge e by 
e* = e/m that allows us to keep all scattering events 
inside the first Brillouin zone. Finally, the distance be- 
tween stripes can be estimated as a y /2 = kp/itn. 

Double strip. Now we consider a double strip, consist- 
ing of two coupled strips for spinless fermions, see Fig. [3] 
This system is equivalent to a single strip but for spin- 
ful fermions. Below we focus on the double strip but we 
note that one can identify the upper (lower) strip with 
spin up (down) state labeled by r\ — I (77 = — 1). The 
chemical potentials /j r; are opposite for the two strips, 
/ii = — \i\ = \[2t x , and are chosen such that the system 
is at half-filling. For the spinful strip the role of is 
played by the Zeeman term, fi n = ngfiBB, arising from 
the magnetic field B along z. Here, g is the g-factor, 
and hb is the Bohr magneton. The inter-strip hopping 
amplitude t z is also accompanied by the phase <j) z arising 
from a uniform magnetic field B2 applied along y, 

H = Y~* / P ""^V t r - ffil 

1L z — / j L z v °n,m,dr,?7 ( -'n,m,CT,?j- \yj 

n 7 m,cr,ri 

The amplitude of B2 is chosen so that (j> z = 
(e/hc)Bia x a z = n. This amounts to apply a total field 
B tot = B + B2 in the yz-plane. Moreover, the same 
H z is generated in the spinful case by a £>2-fiold applied 
along y with an amplitude that oscillates in space along 
x with period 2a x , or, alternatively, by Rashba spin orbit 
interaction^ 

Again, we search for wavefunctions in terms of right 




FIG. 3. Double strip for spinless fermions. The intra-strip 
couplings are the same as in Fig. [1] The inter-strip coupling 
along z is described by the hopping amplitude t z . 
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FIG. 4. Spectrum of a double strip obtained by numerical diagonalization of the tight-binding Hamiltonian H2 = H x + H y + H z 
for the same parameters as in Fig. [2] (a) If \t z \ < t y i [t z /t x = 0.01], there are four edge states at any energy within the gap; 
two of them localized at the left and two at the right edge, (b) If t y \ < \t z \ < t yz [t z /t x = 0.05], there is one zero-energy edge 
state (i.e. with fiat dispersion) at each edge, (c) If \t z \ > t y 2 [t z /t x = 1.5], there are no edge states in the gap. 



and left mover fields defined around two Fermi points, 
kpi = tt/4o, x (upper strip) and k F j = 3ir/4a x (lower 
strip). The linearized Hamiltonian density for this ex- 
tended model in terms of the Pauli matrices r\i acting on 
the upper/lower strip subspace is given by [see Eq. ([3])] 

V.2 = H((T2 -> (J2V3) + tzVi T i ■ (?) 

The resulting spectrum is e l,±,v = 
±yJ(Hv F k) 2 + (tyi +pt z ) 2 , with p = ±1. We note 
that the gap vanishes if \t z \ = t y i or \t z \ = t y 2- The 
closing and reopening of a gap often signals a topo- 
logical phase transition. Indeed, imposing vanishing 
boundary conditions at the edges, we find that there 
are two edge states (one at each edge) at zero energy, 
E = 0, if the following topological criterion is satisfied, 
tyi < \t g \ < t V 2, see Fig. @] The wavefunction of the left 
edge state for t z > is given by (with x = na x ) 

^l=o = (f(x),-if*(x),i(-irf(x),(-ir +i r(x)), 

f(x) = e -ikyay/2 e -(k 2 -+ik Fl )x _ C0S ( kyay /2) 

x e -(*i— «*0» + ism{k y a y /2)er {kl +- lkF ^ x . (8) 

The basis (V'i.i, ipi 1, ipi 1, tpi 1) is composed of wavefunc- 
tions ip v ,a defined at the cr-unit lattice site of the 77-strip. 
The smallest wavevectors k^± = \ti ± t z \/hvF determine 
the localization length of the edge state. We note that 
the probability densities ^^(x)! 2 are uniform inside the 
unit cell. 

If \t z \ < tyi, there are two edge states at each edge for 
E inside the gap; see Fig. H) These states, propagating 
in y, have a momentum k y determined by E. This case 
is similar to one strip with spinless particles discussed 
above. We note that the edge states found here are the 
higher-dimensional extensions of the end bound states 
found in one-dimensional nanowire o 26 ^ 50 and ladders 2 ^. 
For \t z \ > t y 2, there are no edge states; see Fig. 0] Fi- 
nally, the most interesting regime here is t y i < \t z \ < t V 2, 
where there is one zero-energy edge state at each edge, 
see Fig. [4] Such states with flat dispersion are expected 
to be strongly affected by interactions. 

Conclusions. We have studied topological regimes of 
strips with modulated hopping amplitudes in the pres- 
ence of magnetic fields. We found topological regimes 
with chiral edge states at filling factors that correspond 



to integer and fractional QHE regimes. We showed that 
double strips sustain topological phases with mid-gap 
edge states with flat dispersion. Optical lattices^ 2 - seem 
to be promising candidates for implementing directly the 
anisotropic tight-binding models considered here. 

This work is supported by the Swiss NSF, NCCR 
Nanoscience, and NCCR QSIT. 



Appendix A: Fourier transformation 

We introduce the momentum k y via Fourier transfor- 
mation, 

r — - r imk v a v / a -1 \ 

where N y is the number of lattice sites in y-direction. By 
analogy, we introduce a momentum k x , 

Cn, ky ,a = (I/V^E 6 ^"^*^' ( A2 ) 

where N x is the number of lattice sites along x. 



Appendix B: Effective Hamiltonian 

Here we derive the spectrum of H = H x + H y in the 
continuum limit following Refs. [26| and |39|. The annihi- 
lation operator in position space ^(x) close to the Fermi 
level is expressed in terms of slowly varying right [R a (x)] 
and left [L a (x)] movers as 

9{x) = Y}R°(?V kpX + L a {x)e-' lkFX ]. (Bl) 
As a consequence, H x results in the kinetic term, 

H l x m = -ihv F J2 I dx [Rt(x)d x R a (x) 

- L\(x)d x L a (x)l (B2) 
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and H y results in a term that couples right and left 
movers, 

H l y m = t yl J dx [L\{x)Ri(x) + h.c] 

+ t y2 [ dx [e- ik y a » R\(x)Li(x) + h.c] . (B3) 



Here, we used the specific choice of the parameters, 
2kFdx = 4> € (0,7r)- It is this term that leads to res- 
onant scattering and opens the gaps at the Fermi points. 
The Fermi velocity vp depends on the Fermi wavevec- 
tor, hvp = 2t x a x sm(kpa x ). The Hamiltonian density 
H corresponding to H Hn = H l J n + H l y m = J dx 
can be rewritten in terms of the Pauli matrices Tj (ct^), 
acting on the right-left mover (lattice) subspace, ^ = 



L±, Ri, Li), leading directly to Eq. ([3]) in the main 

text. 



Appendix C: Wavefunction of left edge state 

The wavefunction of the state localized at the left edge 
of a spinless strip is given by 

e -k 2 x+ik F x _ e -k lX -ik F x j ! (CI) 

where we suppress the normalization factor. Here, we 
introduced the notations e 10 = {E + iJt^ — E 2 )/t y i and 



h = yjt 2 yl - Eyhv F . 
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